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Abstract 

The binary systems that have been detected in gamma rays have proven very useful to study high-energy processes, 
in particular particle acceleration, emission and radiation reprocessing, and the dynamics of the underlying magne- 
tized flows. Binary systems, either detected or potential gamma-ray emitters, can be grouped in difl'erent subclasses 
depending on the nature of the binary components or the origin of the particle acceleration: the interaction of the 
winds of either a pulsar and a massive star or two massive stars; accretion onto a compact object and jet formation; 
and interaction of a relativistic outflow with the external medium. We evaluate the potentialities of an instrument like 
the Cherenkov telescope array (CTA) to study the non-thermal physics of gamma-ray binaries, which requires the 
observation of high-energy phenomena at difl'erent time and spatial scales. We analyze the capability of CTA, under 
difl'erent configurations, to probe the spectral, temporal and spatial behavior of gamma-ray binaries in the context of 
the known or expected physics of these sources. CTA will be able to probe with high spectral, temporal and spatial 
resolution the physical processes behind the gamma-ray emission in binaries, significantly increasing as well the num- 
ber of known sources. This will allow the derivation of information on the particle acceleration and emission sites 
qualitatively better than what is currently available. 

Keywords: Gamma-rays: observations - Binaries: general - Acceleration of particles - Cherenkov astronomy - 
Radiation mechanisms: non-thermal - Telescopes 
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The many spectacular discoveries made in recent 
years by both satellite-borne (AGILE, Fermi) and 
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ground-based gamma-ray telescopes (H.E.S.S., MAGIC 
and VERITAS) have revealed a variety of new sources 
of high-energy particles in the Galaxy. Among these 
sources we can mention star-forming regions, accreting 
black holes and microquasars, early-type stars with very 
strong stellar winds, young isolated pulsars and their 
nebulae and pulsars in binary systems. The physics of 
particle acceleration and interaction in the complex en- 
vironment of such astrophysical systems is extremely 
rich. The detection of very high energy (VHE) gamma 
rays (£" > 100 GeV) by the current imaging atmospheric 
Cherenkov Telescopes (lACT) from the systems PSR 
B1259-63 iili], LS 5039 H, LS I -h61 303 HH and 
HESS J0632-F057 10], as well as the hint of a VHE flare 
in the black hole binary Cygnus X- 1 13] , provides a clear 
evidence of very efficient particle acceleration in bi- 
nary systems containing compact objects (see e.g. 
In addition, one of the brightest Fermi/LAI sources, 
IFGL JlOl 8.6-5856, has been proposed to be a new 
gamma-ray bina ry iM] that could be associated with a 
H.E.S.S. source Uw . Furthermore, there are other bi- 
nary systems from which VHE emission is expected 
from a theoretical point of view |11]. Although they 
have not yet been detected with the current generation 
of Cherenkov telescopes, their emission at HE gamma- 
rays > 100 MeV) has already been reported in some 
cases (e.g. Cygnus X-3, V 407 Cyg and Eta Carinae). It 
is expected that CTA will find new gamma-ray binaries, 
allowing population studies that will have an impact 
on evolutionary models of high-mass binary systems. 
With a few exceptions, most of the gamma-ray bina- 
ries detected, either accreting or non-accreting sources, 
are all within 3 kpc of the Sun, in a volume equal to 
about ~ 10% of the volume of our Galaxy. Assum- 
ing a uniform distribution, although they should follow 
population I stars with more objects in the spiral arms, 
this is consistent with > 50 or so gamma-ray binaries 
in our Galaxy. This number is also dependent on the 
duty cycle of gamma-ray emission: VHE emission in 
HESS J0632-F057, LS I -h61 303, PSR B 1259-63 is 
strongly dependent on orbital phase and in some sources 
the orbital periods can be (very) long. With a ten times 
improvement in sensitivity, CTA should be able to probe 
for gamma-ray binaries of comparable luminosities up 
to the Galactic center. CTA can thus be reasonably ex- 
pected to detect a couple of dozen gamma-ray binaries. 
The VHE counterparts of LS 5039, HESS J0632-h057 
and (possibly) IFGL 1018.6-5856 were discovered in 
the H.E.S.S. Galactic Plane survey. The ten times more 
sensitive Galactic Plane survey planned for CTA should 
thus enable many discoveries of such systems, which 
are otherwise very difficult to uncover by X-ray, optical 



or radio surveys. A survey of the central portion of the 
galactic plane is planned for the beginning of CTA op- 
eration (see (v^ ), which will pinpoint new gamma-ray 
binaries candidates. 

The study of known and/or new compact binary sys- 
tems at VHE is of primary importance because their 
complexity allows us to probe several physical pro- 
cesses that are still poorly understood. Some of these 
systems are extremely efficient accelerators that could 
shed new light, and eventually force a revision of, par- 
ticle acceleration theory (see e.g. |[l3l] ). The parti- 
cle injection and radiation emission mechanisms in bi- 
nary systems vary periodically due to an eccentric orbit 
and/or interaction geometry changes. This may provide 
information on the location of the high energy particles, 
on the energy mechanism(s) powering relativistic out- 
flows, on the nature of the accelerated particles, and on 
the physical conditions of the surrounding environment. 
The presence of strong photon fields allows the study 
of photon-photon absorption and electromagnetic cas- 
cades. All these processes occur on timescales < 1000 
s, a proper study of which would require at least a 5(T 
(standard deviations) detection for ~ one hour exposure 
times. 

The interaction of binary systems with the Interstel- 
lar medium (ISM) could also be powering a new class 
of TeV sources, which could be resolved/detected with 
enough resolution/sensitivity. For a deep study of the 
processes taking place in compact binary systems we 
need to go beyond the present lACT's capabilities. Be- 
low, we report on examples of numerical simulations 
performed to show how the forthcoming CTA observa- 
tory 1 14] could fulfill these objectives. 

The structure of the paper is as follows: in Sect. [2] 
we outline the stellar gamma-ray source classes that are 
idoneous targets for CTA. In Sect. [3] we introduce key 
questions in high-energy astrophysics that CTA can ad- 
dress and the requirements to achieve them. In Sect. IH 
we present the results of some performance tests of 
the capability of CTA to achieve the aims. Finally we 
present a summary in Sect. \5\ 

2. Binary systems with gamma-ray emission 

2.1. Binary systems with young non- accreting pulsars 

PSR B 1259-63 was the first variable galactic source 
of VHE gamma-rays discovered UJ . It has also been 
detected at HE by AGILE d and Fermi/L AT 
The system contains a 09.5 Ve main sequence donor 
(LS 2883) and a 47.7 ms radio pulsar orbiting its com- 
panion every 3.4 years in a very eccentric orbit (see 
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1 1811 and references therein). Particles are accelerated 
in the shock between the relativistic wind of the young 
non-accreting pulsar and the stellar wind of the mas- 
sive companion star These particles, by 
inverse Compton (IC) up- scattering of stellar UV pho- 
tons should produce VHE gamma rays. The strong 
wind- wind interactions may also produce extended syn- 
chrotron radio emission, as recently reported by 

The other binary systems that have been unambigu- 
ously detected at TeV energies, showing gamma-ray 
flux modulations coincident with their orbital periods, 
are LS 5039 with Pprb ^ 3.9 d Q and LS I -h61 303 
with Porb ~ 26.5 d ll23ll . Fermi has also detected emis- 
sion modulated with the orbital period in both systems 
1 24 l25n . Although the nature (black hole or neutron 
star) of the compact object in LS I -f61 303 and LS 5039 



has not yet been determined ||26|, |27|], both systems 
present some similarities with PSR B 1259-63. T hey 
show variable milli- arc second scale radio structure 1281 
11, similar to that found in PSR B 1259-63. VLBA 
images of LS I -h61 303 obtained during a full orbital cy- 
cle show a rotating elongated morphology |[3Q|i , which 
may be consistent with a model based on the interaction 
between the relativistic wind of a young non-accreting 
pulsar and the wind of the stellar companion ||2Q|i . A 
similar behavior has been observed in LS 5039. This 
system was observed with the VLBA during five con- 
secutive days showing an orbital morphological vari- 
ability, displaying one sided and bipolar structures, but 
recovering the same morphology when observing at the 
same orbital phase 113 IH . The broadband emission from 
radio to VHE gamma-rays of the three sources is vari- 
able and periodic, peaking at MeV-GeV energies. 

For LS I -h61 303 and LS 5039, the GeV and 
TeV emission are well anticorrelated. In particular, in 
LS 5039, the GeV emission peaking around the compact 
object superior conjunction/periastron, and the TeV ra- 
diation around inferior conjunction/apastron. The GeV 
and TeV spectra are also roughly anticorrelated, with 
the GeV emission getting harder for lower fluxes, and 
the TeV emission for higher ones (e.g., 11221 13211 ). In 
both sources, the behavior is more or less compatible 
with radiation produced by IC and moderate gamma- 
ray absorption, processes through which the changing 
geometry along the orbit induces a modulation in both 
flux at TeV and GeV, and spectrum at TeV (and GeV if 
cascades were important). Additional efl'ects like vary- 
ing radiative and adiabatic losses could also afl'ect spec- 
tra and fluxes 113 311 . We note that also IC e- pair moder- 
ate cascading can be important l^, 111, [3 • However, 
the relatively low flux around 10 GeV in both LS 5039 
and LS I -h61 303, below the extrapolation of data and 



model predictions, strongly indicates that the emitter, 
although of likely leptonic-hIC nature, should be quite 
complex (and probably located in the periphery of the 
binary system for LS 5039). 

The most recent addition to the selected group of 
gamma-ray binaries emitting up to very high ener- 
gies is HESS J0632-h057 |37]. The source was ini- 
tially detected by the H.E.S.S. experiment 01, but the 
subsequent non detection by VERITAS excluded it as 
a steady gamma-ray emitter |[38ll . The gamma-ray 
variability was confirmed recently by VERITAS and 
MAGIC which reported an increase of gamma-ray flux 
during 2011 February 7-9 (see Refs. ll39ll and ll40ll . re- 
spectively). The increase in TeV flux coincides with 
the time of a large X-ray peak, that could imply that 
the same population of electrons is producing the X- 
ray and TeV emission similarly to that found in 
LS I -h61 303 1 42] and LS 5039 1330. The source shows 
an X-ray periodicity of 321 ± 5 d which has been as- 
sociated to its orbital period ll4lU43ll . No X-ray pulsa- 
tions have been detected so far ll44ll . In the radio band, 
VLBI observations reveal an extended and variable non- 
thermal radio source ll45ll . with its position being com- 
patible with the BOpe star MWC 148. Overall, these 
results are very similar to the multi-wavelength data 
obtained for the other gamma-ray binaries, although 
HESS J0632-h057 displays a spectral energy distribu- 
tion (SED) one order of magnitude fainter. This may 
provide new constraints on the luminosity distribution 
of this subclass of binary systems in our Galaxy. 

A promising candidate to be a new VHE gamma-ray 
source associated to a binary system has emerged re- 
cently. IFGL J1018. 6-5856 was discovered by search- 
ing for periodicities of Fermi sources, and shows in- 
tensity and spectral modulation at GeV energies with 
a 16.6-day period |9]. Given the variability of the pro- 
posed X-ray and radio counterparts, and the spatial co- 
incidence with an 06V((f)) star, these authors proposed 
this source as a new gamma-ray binary. 

2.2. Microquasars 

Microquasars (MQ) are X-ray binaries (XRB) with 
relativistic jets emitting non-thermal radio emission 
through synchrotron radiation. The gravitational energy 
released by accretion feeds the relativistic radio jets and 
powers the non-thermal emission. Two MQs, Cygnus 
X-3 HStIi and Cygnus X-1 lEil, have been recently 



detected at HE gamma-rays and there is evidence of 
VHE gamma-ray emission associated with Cygnus X- 
1 during a flare H but not with Cygnus X-3 |l5o!l. The 
gamma-ray emission detected by Fermi from Cygnus 
X-3 is modulated with the orbital period of the system 
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and is correlated with the radio emission, which is asso- 
ciated with the relativistic jets. Other MQs, Hke SS 433 
lEill, GRS 1915-105 m or Sco X-1 IHl, have also 
been observed at HE and VHE but they have not been 
detected yet. They show however synchrotron radio 
emission that probes their capabiHty to accelerate parti- 
cles up to relativistic energies, and their transient nature 
makes necessary to look at the sources at the right time 
to detect their high-energy radiation . The presence 
of strong photon fields provided by the companion star 
may off'er a good scenario for the production of gamma- 
rays through IC processes |[55l] . Furthermore, proton- 
proton interactions could also take place in the collision 
between relativistic protons and the stellar wind (s^ . 
However, photon-photon absorption may also take place 
in the inner regions of the system, and the VHE photons 
may be strongly attenuated |[57[[34l [58ll. 

It is noteworthy that, potentially, the non-thermal 
processes leading to radio and high-energy radiation 
in high-mass microquasars can be difficult to distin- 
guish from those in young non-accreting pulsars. As 
pointed out, e.g. in and the complex and 
variable/periodic radio structures found in gamma-ray 
binaries, some with a compact object of unknown na- 
ture, may be also produced by jet disruption or gamma- 
ray absorption and secondary radio emission, respec- 
tively. In addition, the processes underlying or affect- 
ing the high-energy radiation, for instance in LS 5039 
and LS I -h61 303 (as described in Sect. 12. lb , would also 
take place in high-mass microquasars. Therefore, de- 
spite some diff'erences in the non-thermal radiation are 
expected, they are likely to be more of quantitative than 
of qualitative nature, and only high-quality data may be 
able to allow the identification of the underlying engine. 

2.3. Collision of the outflow with the interstellar 
medium 

The termination regions of MQ jets can also gener- 
ate non-thermal emission. The detection of VHE emis- 
sion from those regions would represent a new type 
of gamma-ray emitter. In a scenario similar to that 
found in the Sedov-expansion phase of supernova rem- 
nants (SNRs), MQ jets propagating into the medium 
are eventually decelerated, developing forward/reverse 
shocks into the ISM/jet-ejection. In these shocks 
particles can undergo efficient acceleration and could 
produce gamma-ray emission through IC, relativistic 
Bremsstrahlung and proton-proton collision processes. 
The dynamics of these interactions (see e.g. 116 111 ) and 
the possibility of the production of gamma-ray emis- 
sion ( |[62|] . |[63|] ) have been recently addressed. The 
expected fluxes may depend on the MQ jet power, the 



source age, the external medium particle density profile, 
and the duty cycle if transient, and could reach values 
^E>iOGeV ~ few xlO"^^ erg cm"^ s"^ which is a few 
times below the sensitivity level of the current lACTs. 
These sources may be therefore good targets for the im- 
proved capabilities of CTA. 

Large-scale interactions are also expected in the case 
of binaries hosting a young pulsar, when mixed stellar 
and pulsar winds interact with the environment. The 
expected radiative outcome of this interaction is similar 
to that of microquasar jets 11641] . 



2.4. Colliding wind binaries 

Hot stars can generate strong winds and form collid- 
ing wind binary systems (CWB). Shocks are expected 
to form in massive star binaries, in the region where the 
winds from both stars collide. Non- thermal synchrotron 
emission from the colliding wind region in one source 
has been detected ll65l] . which indicates the presence 
of highly relativistic electrons (see also |'6^). These 
systems may also be embedded in dense photon fields 
where IC losses would be unavoidable, making CWBs 
potential high-energy emitters Jivl [68|] . 

An extreme example is the Eta Carinae system |[69l] . 
Gamma-ray emission has been theoretically predicted 
from this source (see e.g. and the emission has 

been tentatively confirmed recently by the Fermi/LAI 
ItIItS] Sind agile H instruments. The predominant 
GeV emission of Eta Carinae, shown in the top panel 
of Fig. [T] seems to agree with what is expected from 
IC and/or neutral pion decay processes in such type of 
system. At VHE, Eta Carinae has not been detected so 
far 11751] . The reported HE flux levels and the spectrum 
make however this source a good target for CTA, since it 
will provide a significantly improved sensitivity at ener- 
gies in the range 30 to 100 GeV as compared to present 
lACTs. MAGIC observations of WR 146 and WR 147 
produced the first bounds on the high-energy emission 
from Wolf-Rayet binary systems 



2.5. Cataclysmic binaries 

Regarding cataclysmic binary systems, variable high- 
energy gamma-ray emission has been recently detected 
from an optical nova of the symbiotic star V407 Cygni 
[77.1 . According to the source SED shown in the bot- 
tom panel of Fig. [T] the extrapolated VHE gamma-ray 
fluxes are too low for current lACTs sensitivities, and 
even with its improved capabilities, CTA may require 
of relatively long exposure times to detect it. 
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3. Key questions and CTA requirements 
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Figure 1 : Top: Spectral energy distributions (SED) of two FermifLAI 
sources at the position of Eta Carinae: 2FGL J 1045. 0-5941 (red 
points) and PSR J1048-5832 (green points). The best-fit model for 
the average spectrum is shown as a red dashed line (1 73]). Bottom: 
SED of V407 Cyg measured over the period from 10 March to 29 
March 2010. 



3.1. Relativistic outflows in binary systems with young 
non-accreting pulsars 

The rapid rotation of the magneto sphere of 
PSR B 1259-63 generates a highly relativistic wind that 
spins down the young neutron star. The properties of 
such winds are a long-standing issue of astrophysics. 
Constraints on pulsar winds have been inferred from 
observations of their interaction with the surrounding 
medium: the supernova ejecta in pulsar wind nebulae 
(PWN) or the ISM in isolated, fast moving pulsars. In a 
binary system, the pulsar wind interacts with the stellar 
wind of its companion in a periodic fashion. The binary 
nature determines the matter and radiation densities 
to be sampled by the pulsar during its orbital motion, 
and allows the production of relativistic outflows in 
the asymmetric postshock region of the pulsar wind. 
Phase-resolved spectra of binaries powered by a young 
pulsar ofl'er novel tests of pulsar wind models on spatial 
scales ~ 10"^- 10^ times smaller than in isolated objects. 

The study of pulsar winds requires phase-resolved 
light curves and spectra from binary pulsars. The sys- 
tems LS 5039 and LS I -h61 303 have been suggested 
to be powered by a young pulsar ||2Q|i . The orbital peri- 
ods of these systems range from a few days (LS 5039) 
to years (PSR B 1259-63). In LS 5039 the spectrum 
changes drastically between orbital phases of high flux 
as compared to phases of low flux. PSR B 1259-63 and 
LS I -h61 303, on the other hand, are only detected dur- 
ing high flux phases. Information on the low flux phases 
is therefore essential for modeling purposes, and this 
requires the increased sensitivity of CTA. An extended 
spectral range is also crucial to study the connection to 
the MeV-GeV emission (due e.g. to IC emission from 
cascade electrons) and to constrain the high-energy cut- 
ofl's linked to particle acceleration. 

Gamma-ray binaries, in particular LS 5039, present a 
very high acceleration efficiency OH. CTA good 
sensitivity in the range 10-100 TeV will allow us to 
probe the acceleration mechanism, the strength of the 
magnetic field, intrinsic orbital/short timescale variabil- 
ity of the accelerator/emitter, and the emitting particle 
nature. This will be possible since the opacity and or- 
bital dependence of photon-photon absorption, and the 
angular IC dependence, sources of extrinsic variability, 
are smaller in this energy range. In addition, the spec- 
trum of the emission is expected to be different depend- 
ing on whether the emitting particles are protons or elec- 
trons. 

Finally, several studies have tried to synthesize the 
population of high-mass X-ray binaries from assump- 
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tions on the progenitor distributions and their evolu- 
tionary paths. Some include predictions on the num- 
ber of binaries in the spin ning do wn pulsar -hmassive star 



stage as a byproduct ||78|,|79|,|80|]. They are quite uncer- 



tain, with numbers ranging from 100 to 1000 depending 
on the details of binary evolution and on the adopted 
lifetime of this phase. Constraining the population of 
PSR B1259-63-like objects is important as this impacts 
HMXB populations but, ultimately, also the population 
of double neutron stars that are targeted by Virgo/LIGO. 

3.2. The accretion/ejection link in microquasars 

The possible detection of Cygnus X-1 by MAGIC 
may indicate that MQs are able to emit VHE gamma- 
rays. These systems, like active galactic nuclei (AGN) 
and gamma-ray bursts (GRBs), are powered by accre- 
tion. How and where the energy is released are still 
poorly understood. Part of the material is ejected at 
relativistic speeds and this is invariably associated with 
non-thermal emission. Importantly, emission from both 
accretion and ejection processes is detected in MQs. 
Their interplay can be studied when the mass accretion 
rate changes during outbursts. The MQ environment is 
better constrained and the timescales are more accessi- 
ble than in GRBs and AGNs. Emission beyond a few 
MeV in some MQs found by CGRO/COMPTEL is es- 
sentially terra incognita, probably because of variability 
and lack of sensitivity. Exploring the VHE domain and 
connecting it to the known anti-correlations between 
thermal X-ray emission (from accreting material) and 
non-thermal radio emission (from the jet) should bring 
new insights into the accretion/ejection processes. De- 
termining the time variability of the TeV to X-ray spec- 
trum will help to know whether the jets have electron- 
positron or hadronic components, and will help as a con- 
sequence to understand the physics of the jet produc- 
tion. Interestingly, powerful radio flares are likely asso- 
ciated with the production of gamma rays ||47l [SlJ, |82|] . 
As these radio flares occur in a specific area of the hard- 
ness intensity diagram, relatively precise predictions of 
a suitable trigger for CTA can be made. 

Achieving this goal requires the monitoring of se- 
lected sources or target-of-opportunity (ToO) quick re- 
sponse with enough sensitivity to detect flaring episodes 
on ~ hour timescales. The discrete ejections observed 
in X-rays, infrared and radio in GRS 1915-1-105, for 
instance, occur on such timescales, and these could 
also be detected in gamma rays. The full CTA ar- 
ray should detect a 10-mCrab flare in one hour, cor- 
responding to a luminosity above 100 GeV of about 
10"^LEdd for a 10 M© black hole at the Galactic cen- 
ter distance. A significant advantage of CTA compared 



to current lACTs is subarray observations of specific 
sources. This will allow long-term monitoring that is 
currently too time-consuming for presently operating 
lACTs. A 0.1 -Crab flare would be detected in one hour 
by a HESS -like subarray and could trigger follow-up 
observations of the whole array. Relating this emis- 
sion to the known disk/jet/corona phenomenology will 
require on the one hand multi- wavelength campaigns, 
and on the other hand to have the possibility of a quick 
CTA response to ToO's from FermifLAI, X-ray satel- 
lites or radio interferometers like LOFAR and the SKA 
pathfinders MeerKAT, ASKAP and ATA. We note that 
ASKAP and LOFAR are starting all- sky transient mon- 
itoring in 2012. 



3.3. Collision of the outflow with the interstellar 
medium 

MQs are capable to transfer large amounts (~ 10"^^ - 
Iq51 gj,g^ q£ kinetic energy and momentum to the 
surroundings through discrete or continuous ejections. 
Significant non-thermal gamma-ray emission in the 
jet/medium shocked regions could be produced, al- 
though dense ambient conditions are required. ISM 
particle densities can vary from ~10"^ cm"^ (hot ISM 
regions) to ~ 10"^ cm"^ (inside molecular clouds). 
In the case of typical galactic densities ~ 1 cm"^, 
the predicted VHE fluxes above 100 GeV are a few 
xlO"^^ erg s"^ cm"^, accounting for both relativis- 
tic Bremsstrahlung and p-p interactions. The size of 
the interaction structures would be ~1 pc, i.e. ~1 ar- 
cminute for a source located at 3 kpc. The detection of 
VHE emission from the jet/medium interaction struc- 
tures would reveal a new type of VHE emitter, probing 
the acceleration/radiation from large-scale shocks and 
constraining fundamental properties of galactic jets (e.g. 
energetics, magnetic field, etc.). 

To detect the jet/medium interaction region requires 
to improve the present lACTs sensitivity (-fewxlO"^^ 
TeV s"^ cm"^) about one order of magnitude and to have 
an angular resolution of 1 arc-minute. Additionally, the 
possibility of probing the < 50 GeV energy range would 
be extremely helpful, since this is where most of the 
gamma-ray flux is expected to be produced. There are 
no required special timing/monitoring CTA configura- 
tions, since the emission is expected to be steady and the 
source injection variability is supposed to be smoothed 
out at the large termination regions. 

Similar requirements to those of jet termination re- 
gions are met also when observing the termination of 
outflows produced in pulsar binaries. 
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3.4. Colliding winds of massive stars in binary systems 

VHE emission is theoretically expected from binary 
systems with high-mass loss and high- velocity winds. 
These systems display some of the strongest sustained 
winds among Galactic objects and have the highest 
known mass-loss rate of any stellar type. Colliding 
winds of massive star binary systems are potential VHE 
gamma-ray emitters, via leptonic and/or hadronic pro- 
cesses after acceleration of primary particles in the col- 
lision shocks. 

The detection of VHE emission from colliding winds 
requires an improved sensitivity with respect to current 
lACTs. To further study these systems in case of de- 
tection, phase-resolved light curves and spectra would 
be required. Such spectra could give us a clue to un- 
derstanding the physical processes behind the emission, 
since the non-thermal particle distribution strongly de- 
pends on the shock conditions at each orbital phase. A 
low-energy array is favoured; a cutoff at a level of -100 
GeV due to the modest shock velocity and finite size of 
the acceleration zone is predicted, and CTA should be 
able to operate at a lower energy threshold than present 
lACTs. Finally, the orbital distances between the two 
stars (~ 10^^ - 10^^ cm) make the emission region too 
small to be resolved out by the current designs of CTA 
configurations. Angular resolution is therefore not a re- 
quirement in this case. 

4. Performance tests 

Several of the studies herein presented are meant as 
examples of the capability of CTA in comparison with 
current lACTs; thus we take here several of the al- 
ready detected gamma-ray binaries, for which a well 
known set of assumptions can be established. We 
have used Monte Carlo simulations of the sensitivity 
of the array for several possible configurations in order 
to explore the capability of CTA to study binary sys- 
tems. In particular, we have conducted simulations us- 
ing configurations B, D, E, I, NA and NB, and subarrays 
S4-2-120, and s9-2-120. The southern configura- 
tions are optimized for low energies (B), high energies 
(D), balanced with a focus on low energies (E), and bal- 
anced with a focus on medium and high energies (I). 



The first number indicates the number of telescopes, the second 
indicates the type (1 for small size telescope, 2 for medium and 3 for 
large), and the last number indicates the separation in meters. The 
subarray s9-2-120 is composed of nine medium sized telescopes lo- 
cated at a distance of 120m from each other. 



The northern configurations are optimized for low ener- 
gies (NA), and balanced (NB). The subarrays are HESS- 
like arrays (s4-2-120) and expanded HESS-like arrays 
(s9-2-120). All details can be found in Hinton et al., 
this issue. 

4.1. CTA flux error reduction in known TeV sources 

Accurate estimates of the flux, spectral shape, and 
evolution of known TeV sources are very important for 
constraining the physical parameters of the high-energy 
emitting region. This is even more needed when there 
are several parameters that have to be left free when 
fitting data. To explore the CTA capability to derive 
observables and constrain theoretical models, we simu- 
lated the CTA response on LS 5039. This source might 
not be representative of the class of binaries, but will al- 
low us to compare the improvement from CTA data on 
the present generation of lACTs. We based our simula- 
tions on the results obtained by H.E.S.S. on the source, 
simulating the CTA response under similar conditions 
(above 1 TeV). Since the H.E.S.S. data were taken over 
a long time span and under diff'erent zenith angles, the 
energy threshold was not constant. To make a fair com- 
parison with the H.E.S.S. analysis, we treated our simu- 
lated data the same way as the H.E.S.S. collaboration 
did 10]. The simulated counts and the flux normal- 
ization were extracted above 1 TeV assuming an aver- 
age photon index derived from all data: T = 2.23 for 
dN/dE ~ E~^. Based on the H.E.S.S. results we as- 
sumed a sinusoidal shape of the light curve with a period 
of 3.9 d. For each phase point, and using configuration 
I, we then evaluated the light curve and simulated spec- 
tra of the form described above as seen by CTA, for a 
certain observation time in each phasogram bin. By in- 
tegrating these spectra, we got the flux above 100 GeV 
in each bin. We then extrapolated the obtained flux 
value to the integral flux above 1 TeV and propagated 
the error correspondingly. The result for 70 hours of ex- 
posure time of CTA is shown in Figure [21 The improve- 
ment of observations by CTA is clearly visible. The er- 
ror bars are reduced by a factor of -2-4 with respect to 
the H.E.S.S. data points, assuming the average photon 
index of F =2.23. This can be taken as a direct com- 
parison to the published H.E.S.S. results. However, the 
spectral index of the LS 5039 VHE emission changes as 
a function of orbital phase. A variation of the spectral 
index afl'ects the error estimate. Softening the assumed 
photon index (F ~3), CTA data would improve the er- 
rors by up to a factor of 7, whereas a harder spectrum 
(F ~2) results in a minor error reduction of less than a 
factor of 2. We have also performed similar simulations 
but taking a time binning of 10, 20, 50 and 100 bins per 
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full phase period and using a binning of 14 and 28 min- 
utes. With 10 bins per full orbital period the sinusoidal 
shape can hardly be seen, whereas with 20 bins a sine 
function can be fitted to the data points. Taking 50 bins 
of 28 minutes even substructures can be resolved. To 
obtain the same results and similar error bars like those 
reported by H.E.S.S., CTA would only need ~l/6th of 
its observation time, that is, 50 bins of 14 min each, 
pointing towards the possibility of performing a long- 
term monitoring of the global behavior of the source, 
and accessing the duty cycle of the observed features, if 
any. 

Furthermore we studied the minimum observation 
time for CTA to detect the period of LS 5039 in com- 
parison with the H.E.S.S. one. To do so, we simulated 
CTA observations using a sine function over time that 
reproduces the time structure of the H.E.S.S. flux points. 
From the simulated CTA observations we derived flux 
points for each time bin and used those to construct the 
power spectrum of LS 5039. Whereas H.E.S.S. used 
160 bins of 28 minutes (-70 hours in total) to detect the 
3.9 days period of the system, CTA could detect the pe- 
riod with more than 5 cr with only 160 bins of 3 min (8 
hours in total). This would be a significant reduction of 
observation time for CTA. It has to be kept in mind that 
the significance of the period estimation in the H.E.S.S. 
data is larger than 5 cr (i.e. 8 cr), as all the data available 
at the time were used. 

We studied the modulation of the photon index and 
the flux normalisation with the orbital period for a 
source like LS5039. To compare with the H.E.S.S. 
measurements, we assumed 7 hours of observation time 
for each phase bin and simulated the CTA spectra for 
each phase bin with the spectral parameters obtained by 
H.E.S.S. (photon index and difl'erential flux at 1 TeV). 
By fitting these simulated spectra, we obtained the fit 
parameters with the corresponding error. The results are 
shown in Figure [3] The direct comparison of the errors 
of the H.E.S.S. and CTA measurements shows that ob- 
servations with CTA can reduce the errors on the spec- 
tral parameters by a factor between 2 and 4.5. 

The larger sensitivity of CTA would allow tracking 
the behavior of a source in shorter timescales. In partic- 
ular, it would allow comparing with predictions of the 
spectral evolution of a source such as LS 5039, even 
at the minimum of its TeV flux. As an example, we 
used the spectra in phases 0.2 and 0.3 as derived by 
|[83l] . where electromagnetic cascades were included. In 
Figure |4] we show the results of our simulations: in the 
top panel, the two simulated spectra are plotted, assum- 
ing an observation time of 5 hours. Since the recon- 
struction of the energy spectra in true energy requires 




Figure 2: CTA simulations of LS5039. Phasograms of H.E.S.S. 
(gray) and simulated CTA observations (blue). With CTA observa- 
tions the error on the flux can be reduced by a factor 2-4 above 
ITeV. 
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Figure 3: CTA simulations of LS 5039. Top: Photon index versus 
phase (CTA blue, H.E.S.S. gray). Bottom: Flux normalisation versus 
phase (CTA blue, H.E.S.S. gray). 
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Figure 4: CTA simulations of LS 5039. Top: Intrinsic (line) and 
simulated (data points) energy spectra for phases 0.2 (blue) and 0.3 
(black). The assumed integration time is 5 h. Middle: Excess events 
obtained from the simulations above versus reconstructed energy. Bot- 
tom: Residuals of the excess events distributions from the two distri- 
butions above. In all three panels a very clear difference between the 
two spectra can be seen. 



a complicated unfolding procedure we conservatively 
choose to compare the two spectra on the level of the ex- 
cess events as a function of reconstructed energy. The 
two corresponding distributions are shown in the mid- 
dle panel of Figured The distributions are compared to 
each other by calculating the residuals between the two, 
which is shown in the bottom panel. The probability of 
these spectra to be consistent (i.e. to originate from the 
same original distribution) is «:!%. We conclude that 
CTA would easily distinguish between spectra at differ- 
ent phase bins. 

4.2. Short time scale flux variability 

We have studied the capabilities of CTA to de- 
tect short timescale flux and spectral variabilities from 



gamma-ray binaries. The perfect candidate for such 
studies is LS I -h61 303. A simultaneous multiwave- 
length campaign on this source resulted in the discov- 
ery of correlated X-ray/VHE emission with the MAGIC 
lACT and the X-ray satellites XMM-Newton and Swift 
in the energy range 0.3-10 keV |42] (see however [841, 



8511 ). Additionally, fast X-ray variability on timescales 



of a thousand seconds has been detected ||86l |87|], so 
it is natural to think that such fast variability may be 
present in its VHE emission. The detection of such a 
variability would have strong implications on the loca- 
tion and size of the non-thermal emitter. Several theo- 
retical models put forward until now consider the possi- 
bility of short timescale (-100 s), correlated X-ray/VHE 
variability as the result of leptonic emission in small 
scale instabilities and shocks within, e.g. a relativis- 
tic jet Issll or the interaction of the pulsar wind with 
a clumpy stellar wind |[89ll . The detection of such a fea- 
ture with CTA will provide information on the dynam- 
ical mechanisms underlying the VHE emission, which 
is a crucial step in understanding the apparently erratic 
variability of LS I -h61 303. This knowledge will be 
also valuable in understanding other gamma-ray bina- 
ries with complex stellar wind geometries owing to the 
decretion disk of Be stars, such as PSR B 1259-63 and 
HESS J0632-F057. 

Under the assumption that the X-ray/VHE emission 
correlation found by ll42ll holds at shorter timescales, we 
took the X-ray light curve from a ~ 100 ks long Chandra 
observation ll86ll and computed the expected VHE emis- 
sion from the source. We considered the large statistical 
errors in the correlation by varying the parameters fol- 
lowing normal distributions with a width equivalent to 
the uncertainties for each of the points in the light curve. 

We considered two feasible observation scenarios. In 
the first, the full CTA I array is used to perform a deep, 
short observation of LS I -h61 303. In this case, the ob- 
servation length would be of the order of 4 h, which is 
the typical time span for which a source is visible under 
optimum conditions for a single night. To simulate this 
scenario, we selected 10 ks of the X-ray observation and 
obtained the simulated VHE fluxes for energies above 
100 GeV, taking advantage of the improved low-energy 
threshold of CTA with respect to MAGIC. We note that 
the CTA performance simulation is for an observation 
at low zenith angles, and a 4 h observation is likely to 
cover a large range of zenith angles resulting in an in- 
creased energy threshold. 

CTA is able to clearly detect the fast variability exhib- 
ited by the X-ray light curve at scales down to 1000 s, 
as well as recover the X-ray/VHE correlation, with r > 
0.9, for flux variations of a factor 1.5. 
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Figure 5: CTA simulations of LS 1 +61 303 during 100 ks, a time span similar to that of a 10 night observation campaign with around three hour 
of observation per night. Left: X-ray and VHE light curves for CTA s9-2-120. Measured X-ray flux (top, 0.3-10 keV) and simulated VHE fluxes 
obtained with the subarray s9-2-120 (bottom, E > 100 GeV in units of phcm~^ s~^) for the case of a campaign of ten nights of observation and 
binning time of 2500 s. Right: X-ray/VHE correlation for CTA s9-2-120. X-ray/VHE correlation plot from the ten-night, subarray s9-2-120 
observation plotted in 2500 s bins. The Pearson correlation factor, the chance probability of correlation and the parameters of the fitted correlation 
are shown as insets. 




Figure 6: CTA simulations of LS I -1-6I 303. Left: Sensitivity of CTA E to spectral variations. Simulated spectra and power-law fits for 3 hours 
of observation with the full CTA E array for each of the spectral states. The measured power-law indices and the significance of the variations 
are shown. Right: Sensitivity of CTA s9-2-120 to spectral variations. Simulated spectra and power-law fits for 10 hours of observation with the 
subarray s9-2-120 for each of the spectral states. The measured power-law indices and the significance of the variations are shown. 
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The second observation scenario consists of a longer 
observational campaign of about ten nights with three 
hours of observation per night, which we here simulate 
from the continuous lOOks X-ray light curve. Given 
the foreseeable pressure to obtain observation time with 
CTA, this campaign would not be done with the full ar- 
ray but with a so-called subarray, a subset of the 100 
telescopes that may operate independently. Here we 
chose a subarray configuration similar to an expanded 
H.E.S.S.: nine medium sized telescopes located at a dis- 
tance of 120m from each other, known as s9-2-120. 
This kind of configurations will be readily available af- 
ter the construction of CTA begins, and thus the science 
case described could be achieved even before the whole 
array becomes operational. In Figure [5] (left) we see 
that the longer observation time allows the campaign to 
probe a wider range of X-ray fluxes than the variation of 
a factor 1.5 used for the shorter, full CTA array observa- 
tion discussed above. Even with the reduced sensitivity 
of the subarray as compared to the full array, the corre- 
lation is well recovered (see Fig. [5] right), partly owing 
to the aforementioned wider X-ray flux coverage. 

To further explore the shortest time scales in which 
CTA can resolve a flare, we simulated a 20 hours event 
whose flux variation follows a Gaussian distribution 
and assuming the best-fit spectral shape reported by 
MAGIC for the Cygnus X-1 signal H. 20 hours rep- 
resent the total duration of the flare (i.e. the mean 
of the Gaussian distribution is at 10 hours from the 
start of the observation). With a binning of 5 min for 
each data point, CTA could clearly resolve this assumed 
Gaussian-shaped flare. A 5 min integration would result 
in a detection with a significance of Icr at the assumed 
low state and 25(T in high state, whereas with the sensi- 
tivity of MAGIC it is only possible to detect the peak of 
this flare. This is a clear example of the better sensitiv- 
ity of CTA with respect to the existing lACTs. Although 
the limited duty cycle of Cherenkov telescopes prevents 
them to observe a source for 20 hours in a row, a realis- 
tic exposure time of 5 hours would be enough for CTA 
to resolve parts of the flare in bins of ~10 min. 

As a conclusion from this simulation exercise, we 
see that the full CTA array will be a powerful tool to 
probe into the fast flux variability of gamma-ray bina- 
ries, which could allow the characterization of the dy- 
namical processes taking place in the emitting plasma. 
In addition, we have shown that subarray operations 
must be encouraged since they are able to provide very 
interesting results using only a fraction of the full array. 
They will prove invaluable not only in observational 
campaigns where the long-term variability is important, 
but also to monitor flaring sources such as gamma-ray 



binaries and AGNs. 

4.3. Sensitivity to spectral shape variations 

There have been some hints of VHE spectral vari- 
ability of LS I -h61 303 along the orbit [W], but they 
were not significant enough to be claimed as real. The 
proper characterization of the spectral variability would 
allow the characterization of the radiation mechanism, 
and more generally, as mentioned in Section 14.11 the 
determination of the conditions in the high-energy emit- 
ter. To explore the power of CTA to disentangle diff'er- 
ent spectral indices, we have considered a variation of 
the power-law VHE spectrum of LS I -i-61 303 within 
the statistical error obtained in the MAGIC campaign 
on September 2007 (See 1I42I1 ) and test the time required 
by CTA to significantly detect it. The measured spec- 
tral index was r = 2.7 + 0.3, so we took the two ex- 
treme values (2.4 and 3.0) and obtained the simulated 
CTA spectra. We computed the CTA spectra above the 
energy thresholds of 70 GeV and 100 GeV for the full 
array configuration E and the subarray s9-2-120, re- 
spectively. We found that an observation with the full 
CTA E array of 3 hours for each of the spectral states 
(i.e., one night per state) provides a set of spectra from 
which the photon index variations can be clearly seen. 
An example of such a realization can be seen in Fig.O 
left. The diff'erence between the hard (Fh = 2.4) and 
soft (Fs = 3.0) is detected at a confidence level of more 
than 3 cr in 99% of the realizations, with a mean value of 
(7.8 ± 1.8)(T. Considering the variation with the mean 
spectral shape of Fm = 2.7, the significance of the varia- 
tion with respect to the hard and soft spectra is lower but 
still above 3(T for 85% and 74% of the realizations, re- 
spectively. We note that an increase in observation time 
to 5 hours per spectrum leads to a detection rate above 
95% for all the spectrum pairs (see Figure O. We also 
considered the possibility of longer observations with 
the subarray s9-2-120. In this case, for an observation 
of 10 hours for each of the spectral states, the spectral 
variation between the soft and mean spectra is still sig- 
nificant at a level higher than 3(T for 82% of the simu- 
lation realizations. A representative realization of three 
simulated spectra with the subarray s9-2-120 can be 
seen in Fig. 6, right. 

To further test the CTA spectral capabilities, we have 
used the derived spectrum of Cygnus X-1 during the 
flare to simulate 20 energy spectra with photon indices 
ranging from -2 to -4. We have also simulated difl'er- 
ent exposure times: 5, 15, 30 and 60 min, to study the 
minimum time scale to distinguish between the slope of 
difl'erent spectra. Figure [8] left shows the photon index 
error versus the simulated photon index in the fitting of 
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Figure 7: Results of 10000 realizations of the spectral observation 
of LS I +61 303 with the array E during 5h per spectral state. Top; 
Distribution of recovered spectral indices for soft (Fs), mean (Fm), 
and hard (Fh) spectra. Bottom: Distribution of significances of the 
diff'erence between each of the spectrum pairs. The legend indicates 
the fraction of realizations with a significance larger than 3cr. 



the resulting CTA spectra. CTA would be able to dis- 
tinguish the different spectral slopes in all cases except 
those showing the softest spectra, where the error bars 
are too large to properly distinguish them at a high con- 
fidence level. In that case, the observation of a flare as 
that reported in Cygnus X-1 would require exposures > 
15 min. Should such a kind of flare happen again, the 
minimum timescale for a 5 standard deviations detec- 
tion of a flare within 10% of the reported spectrum from 
Cygnus X-1 is in the range of 2-3 min. However, in or- 
der to have a spectrum determination, one can consider 
a lOcr detection threshold. With this constraint, the min- 
imum timescale accessible is in the range of 8.5-12.5 
min. This is shown in Fig[8]right. The above estimates 
hold, provided that the responses of the array are as sta- 
ble as simulated for 30 minutes exposure and that the 
timescales are probed a priori. 

It is clear that CTA will be a powerful tool for the 
detection of spectral variations in gamma-ray binaries. 
The statistical errors we obtained for these simulations 
are nearly an order of magnitude lower than the ones 
obtained with MAGIC in 2007, thus demonstrating the 
capability of CTA to deliver new and exciting science in 
the following years. 

4.4. Exploring the minimum detectable time delay be- 
tween X-ray and TeV emission in gamma-ray bina- 
ries 

As stated in section 14.21 the MAGIC collaboration 
reported on the existence of a correlation between the 



X-ray emission and VHE gamma-ray emission from 
the gamma-ray binary LS I -h61 303 |42] (see however 
1.84.1 ). The electron cooling times owing to synchrotron, 
IC and adiabatic losses inside the system are expected 
to be around a few thousand seconds (see, e.g. 119 111 ). 
If there is a time delay between the non-thermal emis- 
sion at difl'erent bands (e.g.. X-ray and TeV) larger than 
the electron cooling time scale, the emission from these 
bands would most likely have an origin in difl'erent lo- 
cations in the binary. The detection of such a delay has 
been impossible with the current generation of lACTs, 
so here we present a study of the capability of CTA for 
such a study. Furthermore, since CTA will operate to- 
gether with a new generation of X-ray telescopes, such 
as the Japanese X-ray telescope Astro-H, we have done 
this study considering the capabilities of this new in- 
strument. LS I -h61 303 is the binary system selected 
for this study. 

Based on the short X-ray flares from LS I -1-61 303 
detected in fsij], we have chosen to model a Gaussian 
flare with a width of 1500 s. In the Chandra observa- 
tions, these flares have a baseline count rate of 0.35 s"^ 
and a peak of 0.8 s"^ corresponding to unabsorbed en- 
ergy fluxes of ~ 8 and ~ 19 times 10"^^ ergcm"^ s"\ 
respectively. We used these parameters to simulate the 
light curve detected by Astro-H. Using the correlation 
between X-rays and VHE gamma rays found in [42] and 
the tools for simulating the CTA response, we generate 
the corresponding VHE light curves of the flare as seen 
by CTA array above 65 GeV in configuration I. We have 
used a time binning of 600 s for both the CTA and the 
Astro-H light curves (LC), and studied positive delays 
of the TeV light curve with respect to the X-ray light 
curve in the range to 2000 s in steps of 100 s. We 
show in the top panels of Figure [9] the simulated light 
curves corresponding to 500-s delay (left) and 1000-s 
delay (right). In the middle panels of Figure [9] we show 
the X-ray fluxes as a function of the TeV fluxes in both 
cases. A low correlation coefficient is obtained due to a 
loop structure induced by the delay. When the delay is 
the average correlation coefficient is r=0. 86+0.04. 

To clearly detect such delayed correlations, we have 
used the z-transformed discrete correlation function 
(ZDCF), which determines 68% confidence level inter- 
vals for the correlation coeflftcient for running values of 
the delay (see, e.g. Js^, S]). We show in the bottom 



panels of Figure [9] the ZDCF in the case of 500-s delay 
(left) and 1000-s delay (right). The detection of corre- 
lated signals is clear in both cases. Since the VHE LC is 
generated from the X-ray one, their errors are correlated 
and the scattering of one original (X-ray) light curve af- 
fects the second one (VHE). To correct this problem we 
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Figure 8: CTA simulations of Cygnus X-1. Left: Error on the photon index versus the fitted photon index. The photon index error can be understood 
as a proxy of the finest possible resolution. The lines show the fit of a parabolic function to the data points normalized to the first point. Right: 
Integral sensitivity above 30 GeV as a function of integration time for the spectrum of Cygnus X-1 |7], above 5 and 10 standard deviations threshold 
(see text). The shaded area represents a variation of 10% from the spectrum of Cygnus X-1 during the 2007 flare. 



went through simulations: we started with 100 Astro- 
H LC. With each of this 100 LC we produced 10 other 
X-ray LC by adding a Gaussian noise to the original 
ones, for each simulated delay. Then we simulated the 
corresponding VHE LC. At the end we had, for each 
delay, a sample of 1000 pairs of LC. For each pair of 
light curves we have calculated the ZDCF. To evaluate 
how significant is the measurement of a delay using the 
ZDCF we have fitted Gaussian functions to the max- 
ima of the ZDCF. At the end we have 1000 values of 
the peak for each simulated delay distributed around the 
real delay (see Figure [T0|). The measured delay is then 
calculated as the mean value of the distribution and its 
uncertainty is the standard deviation. In Figure [TT] we 
show the measured delay as a function of the simulated 
delay following the procedure described above. Consid- 
ering all possible uncertainties, delays of -1000 s can be 
significantly detected at a 3 cr confidence level in simul- 
taneous light curves obtained with Astro-H and CTA. 
These results are, to first order, independent of the dura- 
tion of the short flares, as far as they last longer than the 
binning. Overall, these results indicate that CTA will al- 
low us to localize and constrain the X-ray ant TeV emit- 
ting regions of gamma-ray binaries and their properties. 

4.5. Exploring the collision ofmicroquasar jets with the 
interstellar medium 
To simulate the CTA response to the observation of 
MQ jet/ISM interactions, we have used the theoretical 
predictions for a source with a jet power of 10^^ erg s"\ 
a source age of 10^ yr embedded in a medium with par- 
ticle density of 1 cm"^ |[62|] . Gamma-ray spectra with 
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Figure 10: Distribution of the fitted peak of the ZDCF for the simu- 
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gives the measured delay (Ar) and its standard deviation cr^^. 
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Figure 9: CTA simulations of LS I +61 303. Top left:X-ray and VHE light curves of a Gaussian flare with sigma of 1500 s and 500-s delay. Top 
right: same with 1000-s delay. Middle left: X-ray flux vs. VHE flux for 500-s delay. Middle right: same with 1000-s delay. Bottom left: ZDCF for 
500-s delay. Bottom right: same with 1000-s delay. 
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Figure 11: CTA simulations of LS 1 +61 303. Measured time delay 
as a function of the simulated delay. Delays as short as 1000 s can be 
significantly measured. 



Fph ~ 2.45 and Fph ~ 2.85 are derived from current 
theoretical models for the leptonic and hadronic contri- 
bution, respectively Ji^, • We studied the CTA per- 
formance in 50 h of observation time using the simu- 
lation tools for the B and E array configurations. Fig- 
ure [12] shows the obtained SED. The simulated flux is 
at a level ~ 1 % of that of the Crab Nebula, although the 
steepening of the spectrum at high energies would make 
it difficult to detect the sources above a few TeV. 

Regarding the extension of the emission as seen in 
gamma rays, accelerated particles emitting at VHE do 
not have time to propagate to large distances since ra- 
diative cooling is very eff'ective. The emission will be 
mostly confined to the accelerator region itself. The 
emitter size may not largely exceed the width of the 
jet, ~ 1 pc, in the reverse shock region. In the case 
of the bow shock, although it may extend sideways for 
much larger distances, only the region around its apex 
(~ few pc) will eff'ectively accelerate particles up to the 
highest energies. The total angular size of the emission 
from a source located 3 kpc away may then be < few 
arcmin, and CTA would image a point-like source with 
only a marginal extension roughly perpendicular to the 
jet direction. 

To resolve the TeV emission produced in the 
jet/medium interaction regions and disentangle it from 
the putative contribution produced close to the jet base. 



the reverse/bow shocks need to be located at a distance 
> 10^^ cm. Although the precise location of the inter- 
action regions for a particular source can be difficult 
to predict since they depend on the jet power and age 
and the medium density surrounding the system, hot 
spots displaying non-thermal emission at similar dis- 
tances have been found, e.g. in SS 433 ll94ll . Upper 
limits to the TeV emission from jet/medium interac- 
tion regions have been already reported for some mi- 
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croquasars, e.g. Cygnus X-1 m and GRS 1915+105 
1 5211 . However, clear evidences of non-thermal emission 
at the interaction sites in those cases are still lacking 
(see, e.g. SSS). 

Steady gamma-ray fluxes are predicted in jet termi- 
nation regions. A constant flux pedestal level could 
be observed on top of the orbital modulations expected 
from the central parts of the system for a wide range 
of orbital geometries and system inclinations due to IC 
and pair-creation angle dependencies in their respective 
cross-sections, and/or due to orbital variations of the ac- 
cretion rate in both hadronic and leptonic scenarios. We 
note that CTA may not be able to separate the possible 
contribution coming from the reverse shock and the for- 
ward (bow) shocks. The shocked jet and medium ma- 
terial are separated by the contact discontinuity, and al- 
though the extent and relative position of it with respect 
to both shocks is difficult to predict (see e.g. ll98ll ). it 
may be too short for the angular resolution at the level 
of few arcminutes expected for CTA. Possible fast dif- 
fusion of accelerated particles behind the shocks would 
complicate further the situation. 

4.6. Exploring the colliding winds of massive star bi- 
nary systems 

We performed numerical simulations of the response 
of CTA for a CWB like Eta Carinae. We based our 
simulations on the measurements of the energy spec- 
trum of Eta Carinae (see top panel of Fig. [T]) by the 
Fermi/LAI ll73ll and the upper limits derived by the 
H.E.S.S. Collaboration ^9^. The spectrum between 0.1 
and 100 GeV is best fit by a power law with an exponen- 
tial cutoff plus an additional power law at high energies. 
In the TeV range. Eta Carinae has not been detected. 
In Figure [13] (left) we show the Fermi/LAT data points 
and the H.E.S.S. upper limits in gray. From these mea- 
surements, it seems that there must be a cutoff' in the 
spectrum at high energies. For our simulations we as- 
sume exponential cutoff's 3.tE= 100, 150, 200 GeV and 
test how well CTA could detect those. We produced 
simulations at increasing observation times in order to 
study the minimal time required to detect the source 



15 



■D 



I : 






source intrinsic spectrum 


- i 


n " decay 








^^^^ CTA simulated spectrum 


* - ■ ^ ■■ 


H H 1 Crab Nebula 
























: 








\ : 




















: \ 








\ 






1 1 1 1 1 1 1 il 1 1 1 1 1 1 1 1 


1 1 1 1 1 iSi\ 1 1 1 1 1 1 II 



0.01 0.1 1 10 100 

E (TeV) 

Figure 12: CTA simulations of the emission from the MQ jet/ISM 
interaction regions. A jet power of 10^^ erg s~^ , a source age of 10^ yr 
and a medium with a particle density of 1 cm~^ have been used. The 
gamma-ray spectra of IC and relativistic Bremsstrahlung and that of 
p-p interactions through tt^ -decay for the adopted parameter values 
are adapted from [62] and |63], respectively. Only the results from 
the B array configuration are displayed, since they are very similar to 
those obtained with the E array. 



and to get a meaningful spectra with such CTA obser- 
vations. In Figure [13] (left), we show the simulated en- 
ergy spectra with different cutoffs as they would be mea- 
sured by CTA. Simulations for 10 hours of observation 
time are displayed. To detect Eta Carinae, CTA would 
need 2-10 hours of observations, depending on the en- 
ergy cutoff in the spectrum; together with Fermi/LAT 
data, it should be possible to determine the cutoff en- 
ergy using a combined fit. However, it would take a 
longer time to determine the cutoff energy using CTA 
data alone. The minimum observation time needed to 
significantly determine the cutoff energy, i.e. to distin- 
guish between a simple power law and a cutoff power 
law, is established using the likelihood ratio test for the 
two hypothesis. In Figure [13] (right), we show the re- 
sulting significance that a cutoff power law is a better 
fit to the data than a pure power law versus integration 
time for the different energies of the cutoff. For this 
study we simulated 100 spectra for each cutoff energy 
and for different observation times as shown in the plot. 
Taking 3(T as a limit to distinguish between the two dif- 
ferent spectral hypothesis, one can see that 20 hours are 
enough to detect the cutoff only if it is above 150 GeV. 
For a cutoff <150 GeV, 30 to 50 hours are needed. From 
our simulations, we can conclude that CTA observation 



times of > 15 hours are necessary to make meaningful 
physics interpretation and modeling, whereas >20 hours 
are necessary to precisely measure the energy cutoff in 
the spectrum. A proper characterization of the highest 
energy cutoff will give important clues on the acceler- 
ation efficiency of the source, which may be operating 
close to the limit predicted by diffusive shock acceler- 
ation, and on the nature of the radiation mechanism, 
either leptonic (IC) or hadronic (proton-proton interac- 
tions). It is noteworthy that other colliding wind binary 
systems hosting powerful WR and O stars may be also 
powerful non- thermal emitters, as hint ed by hard X-ray 
observations or WR 140 with Suzaku lllOOh . 

5. Summary and conclusion 

The sensitivity of CTA will lead to a very good 
sampling of light curves and spectra on very short 
timescales. It will allow as well long source monitor- 
ing using subarrays, still with a sensitivity 2-3 times 
better than any previous instrument operating at VHE 
energies. In particular, it is noteworthy that CTA will 
reduce by a factor of a few the errors in the determi- 
nation of fluxes and spectral indexes. The high sensi- 
tivity and good angular resolution will allow also for 
imaging of possible extended emission in gamma-ray 
binaries, expected at the termination of the generated 
outflows. The low energy threshold will also permit to 
study the maximum particle energy achievable in mas- 
sive star binaries, trace the effects of electromagnetic 
cascades in the spectra of gamma-ray binaries, or catch 
the most luminous part of the spectrum in some sources. 
Finally, under CTA the population of gamma-ray bina- 
ries (and their different subclasses) may easily grow by 
one order of magnitude, which will imply a strong im- 
provement when looking for patterns and trends, tracing 
the physical mechanisms behind the non-thermal activ- 
ity in these sources. For all this, CTA, either in highly 
sensitive observations of the whole array, or under the 
more suitable for monitoring subarray mode, will be a 
tool to obtain the required phenomenological informa- 
tion for deep and accurate modeling of gamma-ray bi- 
naries. This can mean a qualitative jump in our physical 
knowledge of high-energy phenomena in the Galaxy. 
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when compared to a simple power law) versus observation time. The dashed line at 3cr is to guide the eye. 
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